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Abstract. The Oda's Strong Factorization Conjecture states that a proper bi- 
rational map between smooth toric varieties can be decomposed as a sequence of 
smooth toric blowups followed by a sequence of smooth toric blowdowns. This 
article describes an algorithm that conjecturally constructs such a decomposition. 
Several reductions and simplifications of the algorithm are presented and some 
special cases of the conjecture are proved. 



1. Introduction 

The general strong factorization problem asks if a proper birational map be- 
tween nonsingular varieties (in characteristic zero) can be factored into a sequence 
of blowups with nonsingular centers followed by a sequence of inverses of such maps. 
Oda [5] posed the same problem for toric varieties and toric birational maps. Since 
toric varieties are defined by combinatorial data, the conjecture for toric varieties 
also takes a combinatorial form. 

A nonsingular toric variety is determined by a nonsingular fan and a smooth toric 
blowup corresponds to a smooth star subdivision of the fan. The conjecture then is: 

Conjecture 1.1 (Oda). Given two nonsingular fans Ai and A2 with the same 
support, there exists a third fan A3 that can be reached from both Ai and A2 by 
sequences of smooth star subdivisions. 

As the terminology suggests, there also exists a weak version of the factorization 
problem in which blowups and blowdowns are allowed in any order. This weak 
conjecture is known to hold for toric varieties IU [2] and also for general varieties 
in characteristic zero [Hid]. The strong factorization conjecture is open in all cases 
in dimension 3 or higher. 

In this article we study a simple algorithm that is conjectured to construct the 
strong factorization for toric varieties. The problem is that the algorithm may 
run into an infinite loop and never finish. However, computer experiments and 
proofs of several special cases suggest that the algorithm is always finite and solves 
Conjecture II. li 

To construct a strong factorization between two fans, we start with a weak factor- 
ization that is known to exist. In other words, we assume that we can get from one 
fan to the other by a sequence of smooth star subdivisions and smooth star assem- 
blies (the inverses of star subdivisions). The goal is to "commute" star assemblies 
and star subdivisions to have all assemblies after the subdivisions. The algorithm 
takes one star assembly followed by one star subdivision and replaces the pair with 
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a sequence of star subdivisions followed by a sequence of star assemblies. This step 
is repeated until all star subdivisions precede the star assemblies. 

Let us consider 3-dimensional nonsingular fans. We draw such a fan as a cross 
section, which is a 2-dimensional simplicial complex. We may assume that all max- 
imal cones have dimension 3 and that all star subdivisions have their subdivision 
rays in the middle of 2-dimensional cones. In coordinates, a cone generated by 
vi, V2, V3 is divided by a star subdivision into two cones generated by v\ + ^2,^2,^3 
and v±, vi + t>2, V3; the ray generated by v\ + vi is called the subdivision ray. 

Now given two star subdivisions of one fan, we need to construct a common 
refinement by further star subdividing the two new fans. Figure [1] shows two ways 
of doing this, denoted by A and B. In this figure, the fan we start with consists 
of a single cone generated by 1)1,1)2,113. The two subdivisions have subdivision rays 
generated by v\ + V2 and v\ + ^3. Both factorizations A and B replace the assembly- 
subdivision pair by two subdivisions followed by two assemblies. (We read the 
sequence of maps from left to right, starting from the fan Ai and ending with A2.) 




Figure 1. Algorithms A and B. 

Figure Q] describes two factorization algorithms A and B completely. If we have, 
instead of a single cone as in the picture, a global fan and its two star subdivisions, 
then the subdivisions commute if the subdivision rays do not lie in one cone. If the 
subdivision rays do lie in the same cone, then Figure Q] tells us how to factor the 
diagram. (If the cone containing the two subdivision rays lies in a bigger fan, the 
star subdivisions shown in Figure [1] can clearly be exteded to star subdivisons of the 
bigger fan.) 

Figure [2](a) shows algorithm A applied to factor two star assemblies and one star 
subdivision (the lower edge of the diagram) into four star subdivisions followed by 
five star assemblies (the top of the diagram). 

Figure 0(b) shows the first two steps of applying algorithm B to a sequence of two 
star assemblies followed by one star subdivision. One can see that after replacing 
the original cone with the shaded one in the figure, we are back to the situation we 
started with. It follows that algorithm B applied to the two star assemblies and one 
star subdivision in Figure EJb) will run into a cycle and never finish. However, we 
have not found any such infinite loops in the case of algorithm A. Therefore we can 
state: 
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(a) (b) 

Figure 2. Examples of factorization using (a) algorithm A and (b) 
algorithm B. 



Conjecture 1.2. Algorithm A is always finite. 

When constructing examples of factorizations using algorithm A, the most com- 
plicated ones are similar to the one in Figure EKa) . We take two sequences of star 
subdivisions of a single cone generated by fi>^2>^3- On one side we star subdivide 
at the rays generated by v\ + V2 , 2v\ + vi , . . . , mv\ + V2 , and on the other side we 
subdivide at the rays generated by v\ + v?,, 2v\ + v%, . . ., nv\ + V3. The example 
in Figure 0a) shows the case when m = 2 and n = 1. When both m = n = 10, 
the number of star subdivisions in the diagram will be in the thousands. When 
m = n = 40, the number of star subdivisions needed will be in the hundreds of 
thousands. 

One of the main results we prove is that algorithm A is finite on the diagrams 
with m star assemblies and n star subdivisions as described above. We give a 
precise pattern for the cones appearing in the common refinement. On other types 
of diagrams the factorization algorithm may be shorter, but we can not say anything 
about the regularity or patterns that may occur. As a result we cannot prove 
finiteness in general. 

To study algorithm A, we reduce it to the local case and prove that finiteness of 
the local algorithm implies finiteness of the global one. The local algorithm is more 
algebraic. It can be applied to sequences of symbols instead of drawing pictures, 
and it can also be easily implemented on a computer. 

We will work only with fans in dimension 3. The algorithm, however, also applies 
to fans in dimension greater than 3. In higher dimensions we can again assume that 
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all subdivision rays lie in 2-dimensional cones. Then two star subdivisions commute 
unless their subdivision rays lie in the same 3-dimensional cone. In the latter case we 
can apply the algorithm as in the 3-dimensional case. Everything we say below for 
3-dimensional fans is also true, with minimal modifications, for higher dimensinal 
fans. 

Acknowledgements. The two algorithms discussed in this article are certainly 
not new and have been studied by many people. The second author would especially 
like to thank Dan Abramovich, Kenji Matsuki and Jaroslaw Wlodarczyk for fruit- 
ful discussion regarding these algorithms and their possible extensions to non-toric 
situations. 



2.1. Fans and star subdivisions. We refer to [2l[6] for background material about 
fans and toric varieties. 

We only consider 3-dimensional fans in IR 3 where all maximal cones have dimen- 
sion 3. A nonsingular cone a = {vi,V2,Vs) is generated by a basis v\,V2,V3 of the 
lattice Z 3 C I 3 . A nonsingular fan has all its maximal cones nonsingular. A star 
subdivision of a nonsingular fan is called smooth if the resulting fan is again nonsin- 
gular. The inverse of a smooth star subdivision is called a smooth star assembly. 

There are two types of star subdivisions of 3-dimensional fans - the subdivision 
ray can be in the interior of a 3-dimensional cone, or in a 2-dimensional cone. We can 
always replace the first type of star subdivision by a sequence of star subdivisions 
and assemblies of the second type (see Figure [3]). A smooth star subdivision of a 
cone a = (vi,V2,Vs) of the second type has its subdivision ray generated by Vi + Vj 
for i 7^ j. 



Figure 3. Replacing one star subdivision by a sequence of subdivi- 
sions and assemblies. 

2.2. The global algorithm. Recall from the introduction that we start with a se- 
quence of nonsingular fans, connected by smooth star subdivisions. We may assume 
that all subdivision rays lie in 2-dimensional cones. This property is preserved after 
applying one step of algorithm A, hence we will only consider star subdivisions and 
assemblies of this type. The algorithm terminates if all star subdivisions precede 
star assemblies; in other words, when we have a strong factorization. 

At each step of the algorithm there may be many choices of pairs of a star assembly 
followed by a star subdivision that we wish to commute. To make the algorithm 
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not depend on any choices, we need to fix one ordering of such pairs, for example 
we can insist that always the leftmost pair is commuted. However, finiteness of the 
algorithm or its end result (in case it is finite) does not depend on the chosen order. 

It is also clear that the algorithm is finite if it is finite when applied to a sequence 
of star assemblies followed by a sequence of star subdivisions (or just one star sub- 
division). Thus we may consider a single fan A and two sequences of smooth star 
subdivisions of this fan. If the algorithm is finite, it will produce extensions of these 
two sequences resulting in a common refinement. 

2.3. Localization. To localize the algorithm, we replace a fan by a single cone and 
a star subdivision by a subdivision of the cone together with a choice of a cone in 
the subdivided fan. When drawing pictures of local subdivisions, we indicate the 
chosen cone by shading it. 

Now given two local subdivisions of the same cone, we can use the global algo- 
rithm to construct a local factorization. Figure H] shows an example of such a local 
factorization. Notice that the factorization of the two initial local subdivisions in 
this example is unique: there is a unique choice of cone for each subdivision provided 
by algorithm A. 




Figure 4. Local factorization using algorithm A. 

Figure [5] shows two different local factorizations of one pair of initial local subdi- 
visions. In this case algorithm A provides two choices of local factorizations. When 
factoring a sequence that is longer than two star subdivisions, at each step of ap- 
plying algorithm A we have one or two choices of local factorization. As a result, 
we get in general many local factorizations of one initial sequence. 

A local star subdivision can be represented by a matrix as follows. Consider a cone 
(t>i,i>2,f3) and its local star subdivision resulting in the new cone (t>i,t>i + V2,Va). 
Let M be the 3x3 matrix with columns 1)1,1)2,1)3 and N the matrix with columns 
(vi,vi + v 2 ,v 3 ). Then 

N = ME l2 , 

where E12 is an elementary matrix that differs from the identity matrix by the entry 
1 at position (1,2). In the same way all 6 possible local subdivisions of the cone 
(t>i,i>2,f3) can be represented by elementary matrices E^ for i,j € {1,2, 3}, 2 7^ j. 
Star assemblies are represented by inverses E~- 1 of these elementary matrices. 
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A A 
A 7 X A A 7 X A 




Figure 5. Two local factorizations of the same initial sequence. 

To understand the local factorization algorithm in terms of matrices, consider 
Figure [6] where we have labeled the local star subdivisions by elementary matrices. 
The factorization replaces the bottom of the diagram, which we read from left to right 
as E$i, with the top of the diagram E^E^E^ • We denote this replacement 
as: 

E^E^i E^iE^E^ ■ 
(One can recognize this as an actual equality between products of elementary ma- 
trices, but we use the symbol =>• to indicate the direction in which the replacement 
is done.) 



A 

E 32 , . E 12 

A A 




Figure 6. Local factorization represented by matrices. 

The local factorization algorithm can now be described as follows. We start with 
a sequence of elementary matrices and their inverses. At each step we look for a pair 
E^Em in the sequence, and the algorithm tells us how to commute these matrices, 
possibly inserting new elementary matrices. The algorithm is finished when all 
elementary matrices lie to the left of the inverses. 

One problem with the local factorization algorithm that we haven't discussed is 
that the matrix representation M of a cone {1)1,1)2,1)3) depends on the ordering of 
the generators. However, if we choose one ordering of generators, then we automat- 
ically get an ordering of generators after one local star subdivision and hence the 
elementary matrix that represents this subdivision. Now the question is, if we have 
two sequences of local star subdivisions starting and ending with the same cone, do 
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the orderings induced from both sequences agree at the final cone? The answer is 
"no" in general. Figure [7] shows a local factorization step prescribed by algorithm A, 
where there is no consistent ordering of generators in all cones. The two orderings 
of generators in the top cone are cyclic permutations of each other. We write this 
factorization as: 

Here -R321 is a permutation matrix that represents the cyclic permutation of gener- 
ators. 




Figure 7. Ordering of generators in a factorization diagram. 

Figure [8] lists the commutation rules for the local algorithm A in terms of ele- 
mentary matrices. These rules apply for {i,j,k} = {1,2,3}. Rule (1) tells us that 
there is no local factorization and we have to stop the algorithm. The matrices Eij 
and Eji appearing in rule (1) describe the same global subdivision but with different 
choices of cones. Rule (2) tells us to cancel E^ with its inverse. Rules (3) — (6) can 
be read off from the global diagram of algorithm A in Figure [TJ Each of these rules 
corresponds to one cone in the final refinement. Rule (6) gives us a choice between 
two different factorizations; the two factorizations are shown in Figure [5j Rule (66) 
is the only one where there is no consistent labeling of the generators and we need to 
use the permutation matrix R k ji- Rule (7) shows how to commute the permutation 
matrix with the elementary matrices. 

(1) E^Eji => stop. 

(2) EjEij => 1. 

(3) E^. E kj E kj Er. , 

( 4 ) E ij lE jk =>• EjkE^E^ 1 . 

(5) E^Eki => EkiEkjE^ 1 . 
(6a) E^Eik => E lk E^ 

(6b) E^Eik EkiEjkRkjiE^ 1 E~l . 
(7) RkjiEj^ => E±] )r(m) R kji , r : i ^ j ^ k ^ i. 

Figure 8. Rules for local algorithm A. 

As explained above, we start with a sequence of elementary matrices and their 
inverses. The goal is to apply the commutation rules to get all elementary matrices 
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to the left of the inverses of such matrices. We do not care about the location of the 
permutation matrices Rijk', they can be moved to the right or to the left as desired. 
An example of applying the algorithm is: 



At each step we have underlined the pair to which the rule is applied and the rule 
number is shown on the arrow. Note that at the first step we chose to apply rule (66). 
If at the first step we apply rule (6a), then at the next step we would again have a 
choice between rules (6a) and (66). This gives a total of three different factorizations 
of the initial sequence. 

As in the global algorithm, there is a choice of the order in which we apply these 
rules. Since we want to compare the local and the global algorithms, we have to use 
the same order in both algorithms. For example, we can always apply the rule at 
the leftmost place. In the local case when applying rule (6) there is also a choice 
between (6a) and (66). Below, when talking about different choices in applying the 
local algorithm, we always mean the choice between (6a) and (66); we assume that 
the order of applying the rules has been fixed. 

Conjecture 2.1. The local algorithm A is finite: starting with any sequence of 
elementary matrices and their inverses, the rules can be applied only a finite number 
of times for any choices between rules (6a) and (66) that may occur. 

Note that the factorization algorithm ends whenever we apply rule (1) stop, or if 
there are no more places to apply the rules and we have a local strong factorization. 

The main result we prove in this section is: 
Proposition 2.2. Conjecture \2.1\ implies Conjecture li.il 

Proof. Let S be a finite sequence of global star subdivisions and star assemblies. In 
other words, S is a sequence of fans, each one obtained from the previous one by 
one star subdivision or one star assembly. A local subsequence T of S is a choice of 
a cone in each fan of S such that each cone in the sequence is either equal to the 
previous cone, or is obtained from the previous cone by a local star subdivision or a 
local star assembly. A local subsequence gives a sequence of local star subdivisions 
and star assemblies. 

A global sequence S has many local subsequences. In fact, if we choose any cone 
in the first fan in the sequence S, we can extend this choice to a local subsequence 
T. Similarly, if we choose a cone in one fan in the sequence S, we can extend this 
choice to a local subsequence, going in both directions. 

Now suppose S± and S2 are two global sequences and S2 is obtained from S± by 
applying one step of the factorization algorithm. The two sequences can be put into 
a diagram of the same type as in Figure [2J If T\ is a local subsequence of S\ , we seek 
to extend it to a local subsequence T2 of S2 . By this we mean a choice of cones for 
each fan of S2 that agrees with the choice T\ on the fans that are the same in both 
sequences S\ and 52- If such an extension T2 exists, then as sequences of local star 





1 




ON ODA'S STRONG FACTORIZATION CONJECTURE 



9 



subdivisions and assemblies, T2 is either equal to T\ or T2 is obtained from T\ by one 
step of the local factorization algorithm. The extension of T\ to T2 may not always 
exist. However, suppose that instead of T\ we are given a local subsequence T2 of S2, 
then there is always a unique extension of this subsequence to a local subsequence 
T\ of S\. The reason for this is as follows. When we star subdivide a cone, we have 
a choice of two cones to turn it into a local star subdivision; but if we star-assemble 
a cone, there is no choice at all, and it is automatically a local star assembly. If 
T\ is the unique extension of T2, then as sequences of local star subdivisions and 
assemblies, T2 is either equal to T\ or is obtained from it by one step of the local 
algorithm. 

Suppose we have a global sequence of star subdivisions and assemblies S\ on 
which the algorithm is infinite. Applying the algorithm step-by-step, we construct 
new sequences S2, S3, .... Let us also construct a graph with vertices Si and edges 
going from Si to Si+i, indicating that Si + \ is obtained from Si by one step of the 
algorithm. Next we construct a graph of local subdivisions. The vertices are all 
local subsequences of Tj of Si for i > 1 and there is an edge from Tj to Tj + i if Tj is 
a local subsequence of Si and Tj+i is an extension of this to a local subsequence of 
Si + \. By the discussion above, this graph of local subsequences is a set of trees (a 
forest) with roots the local subsequences of S\. Since there are infinitely many *S^, at 
least one of the trees of local subsequences must contain an infinitely long path. We 
claim that one of such infinitely long paths then corresponds to an infinite number 
of local factorization steps applied to a local subsequence of Si, implying that the 
local factorization algorithm is not finite. 

The problem with the claim above is that in the graph of trees some edges corre- 
spond to the identity transformation. It is conceivable that in all infinite paths all 
edges are eventually identities. Since the trees have finite valence, this means that 
for some N > all edges that have distance at least N from the roots are identities. 
To get a contradiction, we can find in some fan in Sm for large M a cone that is very 
"small" in the sense that it takes more than N local subdivisions to reach this cone 
from any cone in a fan in S\ . We extend this choice of a cone to a local subsequence 
Tm of Sm- By construction, this Tm can be reached from a root by more than N 
non- identity edges, which is a contradiction. □ 

The proof of the previous proposition also applies to algorithm B of the introduc- 
tion and its corresponding local version. Since we know that the global algorithm B 
is not finite, the same must be true for the local algorithm. Figure [9] lists the rules 
for the local algorithm B. The rules are similar to the rules of algorithm A. In rule 
(6) we again have a choice between (6a) and (66). In this algorithm there is no need 
for the permutation matrix. A local factorization corresponding to Figure [2](b) is: 



The middle three matrices in the last sequences are the same as in the original 
sequence, hence this algorithm can be repeated cyclically. 

The local rules for the two algorithms are also valid in dimension n > 3. We need 
to assume that {i,j, k} forms a 3-element subset of {1,2, . . . ,n}, and we need one 
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(1) E^Eji stop. 

(2) EtfEij 1. 

(3) E'^Ekj E kj E~- 



( 4 ) E ij E jk E jkE i;j E ik . 

(5) E^E ki EkjEkiE^ 1 . 
(6a) E^E lk ^E ]k E^Ej^. 
(6b) E7j l Eik =>■ EkjEikEr 1 . 



Figure 9. Rules for local algorithm B. 



additional rule 



where i, j, 



k, I are distinct. 



3. FlNITENESS RESULTS 



We prove finiteness of algorithm A in some cases. The term "algorithm" always 
refers to algorithm A. We denote by S, T, . . . sequences of elementary matrices with 
positive powers. We also assume everywhere that {i,j,k} = {1,2,3}. 

3.1. Directed sequences. We call a sequence S directed toward i (or simply di- 
rected) if it consists of elementary matrices E^ and E^. Thus, a directed sequence 
has the form 



Our main goal in this section is to prove that if S and T are both directed, not 
necessarily toward the same i, then the local algorithm is finite on S~ 1 T. 

Let us first add an extra rule to the algorithm that is useful when factoring 
directed sequences. The rule 



allows us to commute Eu with E^,. We claim that rule (8) "commutes" with the 
factorization algorithm in the following sense. Suppose U is a sequence of elementary 
matrices and their inverses and V is obtained from U by applying the rule (8) once. 
Then the factorizations of U are in one-to-one correspondence with factorizations 
of V and and the latter ones are obtained from the former ones by one application 
of rule (8). In particular, if all factorizations of U are finite, then all factorizations 
of V are also finite. To prove this claim, it suffices to compare factorizations of 
E~pEijEik with factorizations of E~^EikEij for different indices a, (3. We do the 
case E a p = Ejf, only and leave the other cases to the reader. 



S = E^E[ 




ik " 



(8) E^ Eik =>■ EikEij. 



Ejk EijEik =>■ EijEikEj k Eik =>■ EijEikEikE - k 
Ejk EikEij E ik E- k l Eij EikEijEi k E- k . 

Clearly the two factorizations differ by one application of rule (8). 
Using rule (8) we can write a directed sequence as 

b ~ h ijE ik - 
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Note also that rule (8) keeps a directed sequence directed and an undirected sequence 
undirected. 

Proposition 3.1. Let S be a sequence directed toward i and T a sequence directed 
toward j, where i 7^ j. Then S~ 1 T has at most one factorization, which is of the 
form TiS^ 1 , where Si is directed toward i and T± is directed toward j. 



Proof. Let 



Then S 1 T factors as 



S 
T 



S _1 T 



E~ n E- m E p -E q , 

%k 13 31 jk 



rpm rpn 

E p -E\. 

31 jk 



E'kE^Er 



stop 



if m = 

if p = 

if m Q,p ^ 



□ 



Now consider the case where S and T are both directed toward i: 

S 
T 



rpm rpn 
^ij-^iki 



To factor 5 _1 T, we can first use rule (8) and rule (2) to cancel elementary matrices 
with their inverses. Depending on the values of m,n,p,q, this brings us to the 
following 4 cases: 



s -l t 



' Tpa rpb 

EijE ik 

rp—a rp—b 
h ij h ik 

Tp—a rpb 
% h ik 

rp—a rpb 



for some a, b > 0. In the first two cases there is nothing more to do. In the last two 
cases, if we only use rule (6a), we can factor: 



rp — a rpb 
&ij ^ik 

E^E\ y 



rpb rp — a 
rpb rp—a 

EijE ik ■ 



From this we get: 

Proposition 3.2. Let S and T be two sequence directed toward i. Lf we do not use 
rule (66) then S~ 1 T has a unique factorization TS" 1 . 

□ 

Corollary 3.3. When rule (66) is removed from the local algorithm A, then the 
algorithm is finite. 

Proof. It suffices to prove finiteness of S~ 1 T, where S consists of one elementary 
matrix, or more generally, where S is directed. Divide T into directed sequences 
T = T1T2 . . . T n . By previous propositions we know that 

< S ,-1 Ti => UV' 1 , 
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where both U and V are directed. By induction on n the factorization of V T2 . . . T n 
is finite. □ 

To finish proving finiteness of the algorithm on 5 _1 T where S and T are directed, 
the only case remaining is when 

Q—lrn rp—m rpn 

b I — E^ £j ik 
and we are allowed to use the full algorithm. 

3.2. Factorization of E~™"E™ k . We will prove below that all factorizations of 
E^ m E^ k are finite. Since we are allowed to use both rules (6a) and (66), there 
are in general many factorizations, and the number of different factorizations grows 
rapidly with m and n. The table in Figure [10] lists the number of different factor- 
izations (not counting the ones ending in stop) for different values of m = n. These 
numbers were found using a computer. 



m=n 


number of factorizations 


1 


2 


2 


6 


3 


16 


5 


68 


10 


658 


20 


8094 


30 


37,322 


40 


112,610 



Figure 10. Number of different factorizations of E-- m E\ 



A group H(j, k, i) is a sequence of elementary matrices of the form 

H(j,k,i) = (EjkEji) 1711 E^j (EjkEji)™ 2 E^j . . . i /•.',/, /•.';, ]'"' /•.',/,/•.',, 

for some 7ij, m,, I > 0. To have a unique expression for H( j, k, i) as above, we require 
that all mi,rii > 0, except possibly mi and ni, and also that I > 0. The shortest 
group is simply H(j,k,i) = Ej^Eij. A partial group H p (j,k,i) is an initial segment 
in a group H(j,k,i). We sometimes write H p (j,k,i) a p to indicate that the partial 
group ends with letter E a p. 

Theorem 3.4. All factorizations of ET: m E™ k are finite and if 

E7™El^TS~\ 
then either T = E™ k or T has the form 

T = E q ik E ki {H{j, k, i)R jkl ) r H p (j, k, i) 

for some q,r > 0. 

Note that since the algorithm is symmetric, S in the statement of the theorem has 
the same form as T (with indices j and k interchanged). The case T = E™ k occurs 
when we commute E~- m with Ef k using rule (6a). The other form of T occurs when 
we apply rule (66) at least once. 
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Let us say that a sequence T has the form o if it is as in the theorem: 

T = E q ik E ki (H{j, k, i)R jki ) r H p (j, k, i). (o) 

Given such a T, we write T a @ to indicate that the last symbol of T is E a p. 
We start with an auxiliary lemma. 

Lemma 3.5. Consider sequences 

(a) T kj E^(H(j, k, i)R jki ) r H p (j, k, i), 

(b) T jk E^R jki (H{j, k, i)R jk i) r H p (j, k, i), 

where T is of the form o. The algorithm is finite on both sequences and produces 
factorizations TiS -1 , where T\ again has the form o. 

Proof. Note that both sequences have a single inverse elementary matrix in them. 
We prove both parts of the lemma simultaneously by induction on the number of 
elementary matrices to the right of the inverse. 
Consider first the sequence (a): 

T kj E^(H(j, k, i)R jki ) r H p (j, k, i) = T k j E^ 1 (Ej k Eji) mi E^ 1 



If mi > 0, the factorization stops with E~^ Eji, so assume mi = 0. If ni > 0, we 
apply one step of the algorithm: 

rp 771 — 1 TP n l T 1 TP TP~^ TP n i~^ 

J-kj&ij E kj ...=> //.•,'•./.•.,'-,, 

We can combine T k jE k j into one Tj£- that again has the form o and we are back 
to the case of (a), but with one less elementary matrix to the right of E^ . If also 
ni = 0, then we apply the algorithm: 

TkjE^ 1 Ej k EijRj ki ... =>• T k jEj k E ik l E^ 1 EijRj ki . . . 

T kj E jkE^ k R jki .... 

We combine T k jEj k into T'- k , and this brings us inductively to the case (b). There 
is also the possibility that in the sequence (a) the group H(J, k, i) occurring is the 
last partial group, and in that group either the last symbol or both Ej k Eij are 
missing. In both cases the algorithm terminates and the form of T\ can be read off 
from the formulas above. 

Now consider sequence (b): 

TjkE^RjkiiHij, k, i)R 3ki ) r H p (j, k, i) => T 3k R 3ki E^(E 3k E T ^E n k ) .... (*) 
When mi > 0, we apply three steps of the algorithm to get: 

rp p TP — If TP TP pill . rp 771 7? ( 771 77! \iril- 1 Z7>ni 

J-jk^jki^ji K^jjk^ji) J^ k j ■■■ => -Ljk^ijTijkiK-^jk-^ji) & k j 

Note that Tj k Eij is of the form o and it ends with a complete group. Thus: 

Tj k EijRj ki = E? k Eki(H(j,k,i)Rjki) r , 
(EjkEji)™- 1 ^ ... = (H(j, k, i)R jki ) s H p (j, k, i). 

Concatenating these sequences gives the T\ as stated in the lemma. 
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When mi = and n\ > 0, we apply the algorithm to the sequence (*) as follows: 
TjkRjkiEjl . . . TjkRjkiEkjEkiEj^E^' 1 . . . 

Combining Tj k EjiEjk = Tj fc , we are inductively back to sequence (6). 
Finally, when mi = ni = 0, we factor the sequence (*) as 

TjkRjkiEj^ EjkEijRjki . . . TjkRjkiEkjEifcRkijE^ E— E^jRjki . . . 

=$> TjfcRjkiEkjEikRkijE^ Rjki • ■ ■ 

This brings us by induction to the sequence (a). We should also consider the case 
where either Eij or both Ej k Eij are missing from the final partial group. Both these 
cases are easy to deal with and lead to the required form of T\ . □ 

Proof of the theorem. To factor E^ m E™ k , we use induction on m. If all factorizations 
of E^ 171 ^ E™ k have the claimed form TS"" 1 , it suffices to prove that all factorizations 

of E~- l T have the same form. The base case m = is trivial. 

v 

When T = E™ k , we can either commute E^ with T using rule (6a), or we can 
commute the first p steps, then apply rule (66): 

E? k E ki E jk R kji E kj Er. E\ k 
=> Ef k E k j Ej k Rkj i E kj E tk E j k E ji 1 

Ef k E k iE jk Rkji {EjkEjjf E^ E jk E ji 1 
=> El k E ki E jk {EjiE jk ) q R kji E kj E jk q Er. . 

Note that Ejk(EjiEjk) q = (Ej k Eji) q Ejk is a partial group H p (j, k, i), thus the result 

T' = Ef k E kl H p (j,k,i) 



is as required. 

assu 

E V lT = B^Kk^iiHU, k, i)R, ki ) r H p (j, k, i 



Now let us assume that T is of the form o and consider factorizations of E--T: 



From the above we know how to factor E^E^. First suppose that the factorization 
is E™ k E~^ . Then we continue with the algorithm: 

ElE^E ki {H{j, k, i)R jkl ) r H p (j, k, i) => ES^E H E k jB^{H(J, k, i)R jkl ) r H p (j, k, i). 

Since the initial segment E™ k EkiEkj = T' k -, where T' is of the form o, we are reduced 
to case (a) of the lemma. 
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Next suppose that we do not commute E^- with all of E™ k : 
> E ik Eki E i k(EjjEjk) 9 Rkj i E^ 1 EJ q ET. 1 E ki H(j, k,i)... 
=> Ef k EkiEjk(EjiEjk) q RkjiE^Ej£ EkiEjj 1 H(j, k,i) . . . 

=> Ef k EkiEjk(EjiEjk) 9 Rkji E k j E k i ( E~ { 1 Ej k l ) q E^ 1 H(j, k,i)... 

= Ef k E ki (E jk E 3i yE 3k R kjl E k lE kl E^(E-^ (**) 

At the next step we cancel pairs Ej k Eji with pairs EJ^EJ^. The number of such 
cancellations depends on q and mi- 
When q < nil — 1, we get: 

(**) => Ef k E k i(Ej k Eji) q Ej k R k jiE k jE k iEj k (Ej k Eji) mi 1 q E k j... 

El k E ki {Ej k Eji) q Ej k R k jiEi k EjiRij k E~^ET k l Ej k {Ej k Eji)^ . . . 



E v ik E ki {E ik E ji ) q E ik Ej i E ki E i ^{Ej k Eji) mx 1 q E k 1 j ... 



Here 



T 1 = Ef k E k i(Ej k Eji) q Ej k EjiE k j = Ef k E k i(Ej k Eji) q+1 E k j 
is of the form o and we are in case (a) of the lemma. 
When q = m\ and n\ > 0, then we get 

(**) E\ k E ki i K Ej k Eji) q Ej k R kji E k ^E ki E~^E'l 1 i . . . 

Ef k E k i(Ej k Eji) q Ej k R k jiE k ^E k iE k jE k iE^ l E k j 1 

Ei k E k i{Ej k Eji) q Ej k R k jiE k iE k ^E k jE ki E-^E k ^ 1 

Ef k E k i(Ej k Eji) q EjkRkjiEjtiEkiEjj 1 E k j 1 . . . 

Ei k E ki (Ej k Eji) q Ej k EijR k jiE ki E-i E k * 1 

In the last sequence we use the fact that R k ji = R 2 j k i and continue: 



Ef k E k i(Ej k Eji) q Ej k EijRj k iRj k iE k iEj i 1 E k j 1 . 
Ef k E ki (EjkEjiYEjkEijRjkiEjkEr^RjkiE^r 1 



This expression falls into case (b) of the lemma. 
When q = mi and n\ = 0, then we get 

(**) => E\ k E k i{Ej k E^ q E^ k R k ^E k jE ki E-^Ej k Eij... 

Ef k Eki(Ej k Eji) q Ej k R k jiE k j E k iE k jEi k R k ijE k ^E i: j Eij . . . 
=> Efk E k j (Ej k Ejj) q Ej k R k jiE k j Ej k R k ijE k l . . . 

Ef k E ki {Ej k Eji) q Ej k EijRj k iRj k iE ik R ki jE k ^Rj ki H(j, k,i) . . . 

=> Ef k Eki(Ej k Eji) q Ej k EijRj k iEf.jE i: i H(j, k,i) ... 

This is the sequence (a) in the lemma. We also have to consider the case where 
the final Em or both Ej k Eij are missing, but these cases are simple and left to the 
reader. 
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When q > mi, we get 

(**) =► ^kiEMEiiE^R^iE^EuiE^E^y^E^ 1 ^ . . . 
When m > 0, this sequence stops with E~ k E k j. When n\ = 0, the sequence 

(**) => fr^(W%^(^ 1 V)'" mi ^- 

=*> E^EkiEjkiEjiEjkyRkjiE^ E ki (E^ Ej k )^~^Ej^ E ih R Uj E^ i E~- . . . 

also stops with EJ k E k j. The cases where either the final E^ or both Ej k Eij are 
missing are left to the reader. □ 

3.3. A global finiteness result. We prove here finiteness of the global algorithm 
A in a special case discussed in the introduction. 

Consider two sequences of global star subdivisions of a single cone (vi, Vj, v k ). The 
subdivivision rays in one sequence are generated by + Vj,2vi + Vj, . . . , mvi + vj, 
and in the other sequence by Vi + v k , 2^ + v^, . . . , ni>j + u& . To prove that the 
algorithm is finite when applied to this sequence of m star assemblies followed by 
n star subdivisions, we follow the notation in the proof of Proposition 12.21 To 
prove finiteness of the global algorithm, it suffices to prove finiteness of the local 
algorithm when applied to all local subsequences of the original global sequence. 
The local subsequences of the m star assemblies - n star subdivisions are: 

TP—m Tpn 

E~ m E q Eu- 
E-- EjJ'Ef, Ei.j, 

ji ij %k Kl ' 

where < p < m and < q < n. Finiteness of the local algorithm when applied to 
the first sequence was proved in the previous subsection. The proof of Thorem 13.41 
also covers the case of the second sequence because T = E\ k E k i is of the form o and 
we proved that for any such T, the factorizations of E~- m T are finite. The case of 
the third sequence follows by symmetry. Only the last case is remaining. Using the 
proof of the theorem, we can factor 

E^E^E? k E ki ^E^TS~\ 

where T is of the form o. Thus it suffices to prove that all factorizations of EJ^T 
are finite. The resulting factorizations may not follow the same pattern as in the 
theorem. 

We apply the algorithm as follows: 

Ej^T = E J ^E^E kl H(j,k,i)... 

> Ef k E^ E^E ki H(j,k,i)... 
=> Ef k E^E ki E^H(j,k,i)... 
=> E q ik E k i( E~i 1 Ej k ) q E^ 1 (Ej k Eji) mi E k j .... 
The part of the sequence that needs to be factored: 

(E^E-^E-^E^E^E^ ... = E^(E^E~lr{E ]k E dl )^El) . . . 
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also appears as a part of the sequence (**) in the proof of Thorem l3.4i We know that 
all factorizations of (**) are finite. This does not directly imply that all factorizations 
of the current sequence are finite. One can, however, repeat the proof of finiteness of 
the algorithm on (**), adjusting the steps where necessary for the current sequence. 

To finish the discussion of the global case, note that in the common refinement of 
the two sequences of global star subdivisions, the number of maximal cones is equal 
to the number of different factorizations of the four types of local sequences above. 
From Figure [10] we know that this number is very large when m and n are large. 
Since each star subdivision of a fan increases its number of maximal cones by one 
or two, it also follows that the number of star subdivisions and star assemblies in 
the strong factorization of the original m-by-ra sequence is very large. 

4. Further directions 

To prove finiteness of the local algorithm A, one needs to find a way to bound the 
complexity of the factorizations, for example their length. The proofs of finiteness 
presented here do not construct such a bound. Instead, we considered very regular 
initial sequences and proved that the factorizations then are also regular. The proof 
of Theorem 13.41 for example does not give a bound on the length of factorizations of 
E~- m E^ k . In computer experiments one can see that the factorizations in these cases 
are in fact rather short. If a factorization has the form TS 1-1 , then the number of 
elementary matrices in T is no more than max{2m + n, m + 2n}. We discuss in this 
section two approaches that may lead to such bounds on complexity and to a proof 
of finiteness. 

4.1. The Cay ley graph. It is well-known that the elementary matrices Eij gener- 
ate the group SL(3, Z). Construct the Cayley graph of SX(3,Z) using these gener- 
ators. This graph has its vertices as the matrices in SL(3, Z) and there is an edge 
from X to Y if Y = XEij for some If we represent cones by matrices, then a 
sequence of local star subdivisions is simply a path in this graph. To make the repre- 
sentation independent of the chosen order of generators for a cone, we should really 
consider the quotient of the Cayley graph by the alternating group A% that acts by 
cyclically permutating columns of matrices. Then all factorization diagrams, such 
as the example in Figure [U can be thought of as being subgraphs of the quotient 
graph. 

Cayley graphs of groups such as SL(3, Z) are studied in combinatorial group the- 
ory. There are many techniques and results known that are similar to our problem. 
For example, in studying the isoperimetric inequalities for the Cayley graph in the 
word metric, one starts with a closed loop and asks how many relations are needed to 
contract this loop to a point. This can be compared to the factorization problem as 
follows: if we start with a partial factorization diagram, which is a loop in the graph, 
we seek to expand this loop by applying the relations. It may be possible to apply 
techniques from combinatorial group theory to get invariants for the factorization 
problem. 

4.2. Factorization along a valuation. Define a valuation as a ray in M 3 generated 
by a vector v with rationally independent coordinates. In the local algorithm the 
valuation ray tells us which cone to choose after a star subdivision - we always 
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choose the cone containing the ray. Using the local algorithm when the choice of 
cones is given by a valuation ray is called factorization along a valuation. 

One can pose the following conjecture: 
Conjecture 4.1. Algorithm A is finite along any valuation 

As in the local case, we also have the following problem: does finiteness of the 
algorithm along any valuation imply finitness of the local or global algorithm? 

The factorization algorithm along a valuation is easier to visualize than the global 
or the local algorithm. If a cone (vi,V2,v 3 ) contains the vector v, write 

v = bivi + b 2 v 2 + b 3 v 3 

for some positive numbers bi,b 2 ,b 3 . We can then represent the cone by the vector 
(61, b 2 ,b 3 ) in l^ 3 - A local star subdivision of this cone along the valuation corresponds 
to subtracting 6j from bj for some i 7^ j, for example 

(bi,b 2 - h,b 3 ) -> (bi,b 2 ,b 3 ) 

is one star subdivision, provided that b 2 > b\. Using this representation of cones as 
points in R 3 , one can also consider embedding a factorization diagram in R 3 . As 
before, to make this independent of the order of generators, the triple (bi,b 2 ,b 3 ) 
should be considered up to cyclic permutation. Because of this cyclic permutation, 
the ambient space M. 3 /A 3 becomes more complicated to work with. In R 3 or in the 
quotient R 3 /^ where one can measure lengths and distances, it may be possible to 
find numerical invariants that bound the complexity of the algorithm. 

We also remark that the local algorithm B does not need the cyclic permutations, 
hence a factorization diagram can be embedded in R 3 . However, it is not hard to 
find a valuation and a diagram such that the factorization algorithm B along the 
valuation is not finite. (The example in Figure Wfo) is actually finite along any 
valuation, but if one takes a symmetric 2-by-2 initial sequence, then there are many 
valuations along which the algorithm is not finite.) 
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